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The wave function of a dilute hard sphere Bose gas at low temperatures is discussed, emp hasizing 
the formation of pairs. The pair distribution function is calculated for two values of \J pa?. 
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The simplest nontrivial quantum mechanical many 
body problem that is precisely definable and at the same 
tim.e subject to mathematical treatment is the dilute 
hard sphere Bose gas. In 1947 Bogoliubov studied this 
problem [l| and obtained its excitation spectrum. His 
method is now called [3| the Bogoliubov transformation. 
In 1957 Lee, Huang and Yang made !3| a detailed study 
of the problem and obtained its ground state energy in 
an asymptotic expansion: 
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where N is the total number of particles, p = N/^ is the 
density of particles, and a is the diameter of the sphere 
= scattering length. 

With incredibly advanced modern technologies this 
system is now subject to experimental study. Many beau- 
tiful results [J| have been obtained, about equation (1) 
and about the excitation spectrum. We want to point 
out in this brief note that perhaps further interesting ex- 
perimental exploration could be made to further clarify 
properties of the ground state. Throughout this note we 
follow the notation of reference [3] . 



WAVE FUNCTION FOR GROUND STATE 

LHY had given a physical picture of the ground state 
^0- It turns out to be a collection of n pairs of particles 
floating in a sea of k = particles. Each pair consists of 
two particles with momenta k and — k. The pair is not 
a bound pair, but a correlated pair with a correlation 
length equals to 
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n = 0, 1, 2..., with n = being dominant, n 
dominant, etc. The average value of n is 
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Thus the probability of a particle being in k 7^ state is 
■:7^ y^pa^. Within one correlation length the number of 
particle is, on the average, 

-pr',^^^»l. (4) 



But most of these would be in the k = state, with only 
a few particle having k 7^ 0: 
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In other words, within one correlation length, there are 
many particles, but only a few with k 7^ 0. The cor- 
relation therefore acts over a sea of particles almost all 
of which have k = 0. That characteristic is, we believe, 
what was vaguely called momentum space ordering by 
London [5| in 1954 for a superfluid system. 

It is interesting to compare the present theory for 
Bosons with the BCS theory for Fermions. In both cases 
there is formation of correlated pairs. But the mecha- 
nisms for the formation are quite different. The impor- 
tant common feature for the two cases is the presence of 
ODLRO 11. 



PAIR DISTRIBUTION FUNCTION 



The pair distribution function 

D{ri2) = p-^{^/j\ri)^/j\r2Mr2Mri)) (5) 

is an important physical quantity measurable for many 
liquid systems. In this formula 4>(r) is the annihilation 
operator in r space. 

The pair distribution function D{r 12) is also related to 
the diagonal elements of the reduced density matrix [61] : 

Diri2) = p~^TTace[^|J{r2Mrl)pN^^irl)^p^r2)], 

where pN = ^0*0 • 

The meaning of D(r) is: given a particle at a point, 
the average number of particles in d^r at a distance r is 
pD{r)d^Y. Thus D{r) ^ 1 as r ^ 00. 

For the dilute Bose system under consideration, D{r) 
has been [7| computed in LHY: 

D{r) = [1 + G(r)]2 + [1 + F{r)Y - 1 - 2/[G(r) + F(r)] 
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where 
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[44] 
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FIG. 1: (color online) Plot of D{r/ro) 



The quantity a was defined in LHY by [38] . All equa- 
tions with [ ] refer to equations in LHY. Converting the 
integrals / rf^k into radial and angular parts yield 



F{r) = ^J'-^hiA), 
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and a == 1 + C^ - ^(1+^^)2- 

We plot [7] D(r) for several values of 
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